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ABSTRACT 

Spectroscopic and photometric observations show that many globular clusters host multiple stellar 
populations, challenging the common paradigm that globular clusters are "simple stellar populations" 
composed of stars of uniform age and chemical composition. The chemical abundances of second- 
generation (SG) stars constrain the sources of gas out of which these stars must have formed, indicating 
that the gas must contain matter processed through the high-temperature CNO cycle. First-generation 
massive Asymptotic Giant Branch (AGB) stars have been proposed as the source of this gas. In a 
previous study, by means of hydrodynamical and N-body simulations, we have shown that the AGB 
ejecta collect in a cooling flow in the cluster core, where the gas reaches high densities, ultimately 
forming a centrally concentrated subsystem of SG stars. In this Letter we show that the high gas 
density can also lead to significant accretion onto a pre-existing seed black hole. We show that gas 
accretion can increase the black hole mass by up to a factor of 100. The details of the gas dynamics 
are important in determining the actual black hole growth. Assuming a near-universal seed black hole 
mass and small cluster-to-cluster variations in the duration of the SG formation phase, the outcome 
of our scenario is one in which the present intermediate-mass black hole (IMBH) mass may have 
only a weak dependence on the current cluster properties. The scenario presented provides a natural 
mechanism for the formation of an IMBH at the cluster center during the SG star-formation phase. 

Subject headings: globular clusters: general 



1. INTRODUCTION 

If globular clusters follow the same u Mbh ~ a " rela- 
tion between central black hole mass and effective ve- 
locity dispersion as is observed in galaxies (Tremaine 
et al. 2002), many clusters should host intermediate- 
mass black holes (IMBHs) with masses in the range 
10 2 < Mbh/M q < 10 4 . Unfortunately, hard observa- 
tional evidence for IMBHs in globular clusters has proved 
elusive. Currently the strongest case is the massive clus- 
ter Gl in M31, where dynamical, X-ray, and radio studies 
are consistent with a mass Mbh ~ 2 x 10 4 (Gebhardt 
et al. 2005, Ulvestad et al. 2007). In our own Galaxy, a 
measurement of an IMBH mass Mbh ~ 4 x 10 4 M q has 
been reported for w Centauri (Noyola et al. 2008), al- 
though this result has been questioned in a recent study 
(Anderson & van der Marel 2010, van der Marel & An- 
derson 2010) suggesting that the IMBH signature be- 
comes much weaker when a more accurate estimate of 
the cluster center is used. Recently Ibata et al. (2009) 
have reported evidence of a density and kinematic cusp 
in the core of M54, a cluster located at the center of the 
Sagittarius dwarf galaxy; their models suggest that the 
cusp could be due to a - 10 4 M Q IMBH. 

For most Galactic globular clusters, kinematic and 
structural data currently provide only upper limits 
on possible IMBH masses (McLaughlin et al. 2006, 
Pasquato et al. 2009, van der Marel & Anderson 2009). 



However, these upper limits are larger than the values of 
Mbh expected from the Mbh — cr relation, and do not 
exclude the possibility that many globular clusters host 
IMBHs with masses in the above range (van der Marel 
& Anderson 2009). 

The presence of an IMBH in a globular cluster has 
important consequences for the cluster's structure, kine- 
matics, internal energetics and long-term dynamical evo- 
lution (see e.g. Baumgardt et al. 2004, Heggie et al. 
2007, Trenti et al. 2007, Miocchi 2007; see also McMil- 
lan 2008 for a review). In addition, Eddington accretion 
onto IMBHs has been proposed as the mechanism pow- 
ering some ultraluminous X-ray sources (e.g. Farrell et 
al. 2009), and IMBHs have long been recognized as po- 
tential rich sources of gravitational waves detectable by 
LISA (Miller 2009). For these reasons, despite the lack 
of strong observational constraints, the formation and 
consequent properties of IMBHs have become a topic of 
considerable theoretical and observational interest. 

One possible IMBH formation mechanism is direct col- 
lapse of a very massive (> 250 Mq) Population III star 
(see e.g. Bond et al. 1984, Madau & Rees 2001 and refer- 
ences therein; but see Abel et al. 2002 for hydrodynami- 
cal simulations suggesting that these stars would form in 
isolation). For non-primordial or more massive IMBHs, 
the leading formation mechanisms center on dynamical 
processes in dense star clusters. These include runaway 
mergers of massive stars in clusters with high central den- 
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sities (Portegies Zwart & McMillan 2002, Giirkan et al. 
2004, Portegies Zwart et al. 2004), and merging of black 
holes in binaries (see O'Leary et al. 2006 and references 
therein). However, numerous potential problems affect- 
ing these models have been pointed out in the recent 
literature, suggesting that the net growth rate and hence 
the final mass of the resulting IMBH may be much lower 
than suggested by earlier dynamical simulations (Yun- 
gelson et al. 2008, Glebbeek et al. 2009). 

Thus, while there appear to be numerous plausible stel- 
lar, binary, and dynamical pathways to the formation of 
relatively "low-mass" (~ 100 - 200M©) IMBHs in clus- 
ters, there is currently no clear consensus on a mechanism 
for the production of high-mass (10 3 - 10 4 M Q ) IMBHs. 

In this Letter we propose a new scenario for the growth 
of IMBHs in early globular clusters. It is a natural conse- 
quence of the physical conditions that may have obtained 
in cluster cores during the formation of the second gener- 
ation (hereafter SG) stars now observed in many globular 
clusters. Specifically, we show that accretion by a pre- 
existing seed black hole of just a small fraction of the gas 
from which the SG stars formed can lead to black holes 
in the high-mass (10 3 - 10 4 M Q ) IMBH range. 

The structure of this Letter is as follows. In Section 
Owe briefly review the observational evidence for multi- 
ple stellar populations in globular clusters, and summa- 
rize the main elements of our model for the formation 
and evolution of multiple populations. In Section [3j we 
present our scenario for IMBH growth. We discuss some 
consequences of this scenario in Section 21 

2. MULTIPLE STELLAR POPULATIONS IN GLOBULAR 
CLUSTERS 

Spectroscopic and photometric observations provide 
strong evidence for multiple stellar generations in globu- 
lar clusters. 

Spectroscopic studies have shown that, even in clus- 
ters whose stars have a common Fe abundance, there are 
substantial star-to-star differences in the abundances of 
lighter elements (see Gratton et al. 2004 and references 
therein, and Carretta et al. 2009a, b for a recent spec- 
troscopic study of 19 Galactic clusters). All the Galac- 
tic globular clusters for which spectroscopic studies have 
been carried out show evidence of multiple stellar popu- 
lations, and in all cases SG stars make up a significant 
fraction (50 - 80%) of the total (D'Antona & Caloi 2008, 
Carretta et al. 2009a,b). 

Several recent photometric studies have provided fur- 
ther support for these results by showing the photomet- 
ric fingerprints of multiple populations (multiple main 
sequences, significant main sequence broadening, or sub- 
giant branch splits) in the color-magnitude diagrams 
both of "normal" clusters in the Galaxy (Piotto et al. 
2007, Marino et al. 2008, Milone et al. 2008, Anderson 
et al. 2009), and of more complex systems such as ui Cen 
(Lee et al. 1999, Pancino et al. 2000, Bedin et al. 2004, 
Villanova et al. 2007), M22 (Marino et al. 2009, Lee et 
al. 2009) and probably Terzan 5 (Ferraro et al. 2009). 

The differences in light-element abundances suggest 
that SG stars formed out of gas containing matter pro- 
cessed through high-temperature CNO cycle reactions in 
first-generation (FG) stars. The main candidates cur- 
rently suggested as possible sources of gas for SG forma- 
tion are massive AGB stars (Ventura et al. 2001) and 



rapidly rotating massive stars (Decressin et al. 2007). 
Although many aspects of these two scenarios are still 
under investigation, in order to form the large mass of 
SG stars suggested by observations, both scenarios re- 
quire that either the IMF of the FG stars was highly 
anomalous, with an unusually large fraction of massive 
stars, or the FG population had a normal IMF but was 
initially at least ten times more massive than is now ob- 
served. 

We have recently carried out hydrodynamical and N- 
body simulations to explore the formation and dynamical 
evolution of multiple population clusters formed out of 
AGB star ejecta (D'Ercole et al. 2008). Our simula- 
tions show that the gas lost by FG AGB stars collects in 
a cooling flow in the cluster core, where most SG stars 
are formed. We have further demonstrated a dynamical 
mechanism whereby, after the centrally concentrated SG 
stars have formed, the cluster's response to early mass 
loss (of intracluster gas and/or supernova ejecta) leads 
to the escape of a substantial fraction of FG stars. Dur- 
ing the subsequent long-term evolution, the SG and FG 
populations mix by two-body relaxation, plausibly lead- 
ing to the multiple-population clusters observed today. 

We now explore the implications for the growth of a 
pre-existing seed black hole of the physical conditions in 
the central regions of the cluster during the SG formation 
phase. 

3. GROWTH OF A SEED BLACK HOLE 

In the hydrodynamical simulations presented in 
D'Ercole et al. (2008), we assumed that the cluster FG 
stars followed a King (1962) density profile 
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out to a truncation radius R t , with stellar masses dis- 
tributed in the range 0.1 < m/rriQ < 100 according to a 
"Kroupa" initial mass function (Kroupaet al. 1993). We 
adopted a concentration c = \og(R t /r c ) = 1.5 and fol- 
lowed the evolution of the AGB ejecta in clusters having 
different values of the total FG mass M and truncation 
radius Rt- 

Figure 1 shows the time evolution of the AGB ejecta 
density and sound speed within 0.1 pc of the cluster cen- 
ter during the SG formation phase, for four clusters span- 
ning broad ranges in initial mass and tidal radius. In all 
cases, if the cluster already hosts a seed black hole with 
mass Mbh ~ 10 2 Mq — the result of one of the "prompt" 
low-mass IMBH formation mechanisms discussed previ- 
ously, operating in the FG population — the gas condi- 
tions shown in Figure 1 will lead to accretion. The Bondi 
accretion rate onto a black hole of mass Mbh is 
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where the density p and the sound speed c s have been 
normalized to the typical values found in our simulations 
and shown in Figure 1. For the physical conditions in our 
four simulations, the Bondi accretion rate is ~ 10 -4 — 
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Fig. 1. — Time evolution of the AGB ejecta density (lower panel) 
and sound speed (upper panel) within 0.1 pc of the cluster cen- 
ter during the SG formation phase for a cluster with FG total 
mass M and tidal radius Rt equal to (10 7 Mq, 200 pc) (solid line), 
(10 7 M ,4Opc) (dashed line), (10 6 Mq, 90 pc) (long-dashed line), 
(10 6 M Q , 18 pc) (dot-dashed line). 



10 3 M©/yr, much larger than the Eddington rate, 
- 3 x 1CT 9 (Mbh/Mq)/^ M Q /yr (Krolik 1999), where 77 
is the radiative efficiency, to be discussed further below. 
The black hole will therefore accrete gas at the Eddington 
rate. 

Under these circumstances, the black hole mass will 
grow on a Salpeter time scale 

T s ~ 3.9 x 10 8 t) LE/Lyv, (3) 

where L E ~ 1.45 X 10 38 M BH /M Q erg/s is the Eddington 
luminosity and a value of the mass per electron [i e ~ 1.15 
has been adopted (see e.g. Krolik 1999). The value of 77 is 
uncertain and depends on the properties of the accretion 
disk and on the spin of the accreting black hole. Here 
we adopt r\ ~ 0.05 — 0.1, a range commonly adopted 
in studies of black hole growth in galactic nuclei. We 
note that values of r\ larger than 0.1, appropriate for 
rapidly rotating black holes, would reduce the growth of 
the IMBH, while much smaller values (i| < 1 ) typical 
of thick accretion disks would significantly increase the 
IMBH growth rate. For accretion at the Eddington rate 
with 77 ~ 0.05 - 0.1, T s ~ 20 - 40 Myr. 

The duration of the SG star formation episode in the 
models of D'Ercole et al. (2008) is AT w 10 2 Myr. This 
time scale is constrained by spectroscopic observations 
showing little variation in the total CNO of FG and 
SG stars, implying that only AGB progenitors with 
M > 4 — 5M© must have contributed gas to the SG 
formation process. Thus, if the seed black hole accretes 
at close to the Eddington rate during this time, its mass 
will grow by a factor of ~ 10 — 100, leading, by the end 
of the SG formation phase, to a final black hole mass of 



Mbh ~ 10 3 - 10 4 M ( 



4. DISCUSSION 

Several aspects of the scenario just presented merit fur- 
ther discussion and exploration. 

We consider first the effect on the cooling flow of the 
energy radiated by the accreting black hole. To address 
this issue we follow the analysis of black hole outflow 
presented by King (2003, 2005), in which it is shown 
that the wind produced by the black hole sweeps up the 
surrounding gas into an expanding shell. Because of ef- 
fective energy losses due to inverse Compton scattering, 
the shell expands in the momentum-conserving regime. 
Under these conditions (see King 2003, 2005), in order 
for the momentum flux from a black hole accreting at the 
Eddington rate to overcome gravity and expel the shell 
from the cluster, the black hole mass must be larger than 

_ fgKO 4 

MBH,crit — — J^T- iU 

Here, k is the electron scattering opacity, f g is the frac- 
tion of the total cluster mass in the form of gas, and it has 
been assumed that the gas is embedded in an isothermal 
system with velocity dispersion a. 

Assuming a > 20 km/s (the smallest value for which 
the AGB ejecta can be retained) and f g ~ 0.05 (D'Ercole 
et al. 2008), we find M BHxrit > 1O 4 M , much larger 
than the assumed mass of our initial seed black hole. 
Thus feedback from the accreting black hole is not suf- 
ficient to alter significantly the global dynamics of the 
cooling flow or the SG star formation process. 

On the other hand, the effects of feedback may signifi- 
cantly alter the local gas dynamics in the vicinity of the 
black hole. For example, recent detailed 2-D hydrody- 
namical models of accretion onto a 100 M© seed black 
hole in a dense protogalactic cloud (Milosavljevic et al. 
2009) show intermittent accretion, reducing the net ac- 
cretion rate to about 1 /3 of Eddington. Such a reduction 
in the mean accretion rate in our case would reduce the 
growth of the seed black hole to a factor ~ 2 — 5 for the 
range of r\ values adopted in Sj3l 

We note here that during the accretion phase the 
black hole will be very luminous and might be observed 
as a (possibly intermittent) ultraluminous X-ray source 
(ULX). We caution however that, before linking this sce- 
nario to the ULXs observed in nearby young massive 
clusters (see e.g. Portegies Zwart et al. 2010), one must 
first verify whether these clusters (1) meet the conditions 
to form a seed black hole, (2) are sufficiently massive and 
(3) have the structural properties required for SG star 
formation. 

A second important issue is the possible relation be- 
tween the IMBH mass and the current structural prop- 
erties of the parent cluster. 

Any correlation between the mass of the seed black 
hole and the early properties of the cluster should be 
preserved during the Eddington growth phase, provided 
that the duration AT of the SG star formation episode 
is approximately independent of the cluster environment. 
The mass of the seed black hole depends on its forma- 
tion mechanism. If it is the result of stellar evolution 
in a massive Pop. Ill progenitor, a roughly universal 
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mass seems the most likely result (see e.g. Madau & 
Rees 2001). Early dynamical simulations of the runaway 
merger scenario suggested that the mass of the runaway 
should scale with the mass of the cluster (or the cluster 
core) in which it forms (Giirkan et al. 2004, Portegies 
Zwart et al. 2004). However, as mentioned in JEJ strong 
stellar winds may well limit the actual mass attained, 
possibly to as little as a few hundred solar masses, again 
largely independent of the initial cluster mass. 

For the black hole growth phase to occur, the cluster 
must initially have been massive and/or compact enough 
to retain the AGB ejecta that subsequently flowed into 
the cluster core and formed the SG population. As dis- 
cussed in D'Ercole et al. (2008), multiple population 
clusters must have undergone an early phase of strong 
FG mass loss and structural evolution; it is therefore 



Abel, T., Bryan, G. L., Norman, M. L., 2002, Science, 295, 93 
Anderson, J., Piotto, G., King, I.R., Bedin, L.R., Guhathakurta, 

P., 2009, ApJ, 697, L58 
Anderson, J., van der Marel, R., 2010, ApJ, 710, 1032 
Baumgardt H., Makino, J., Ebisuzaki, T., 2004, ApJ, 613, 1143 
Bedin, L.R., Piotto, G., Anderson, J., Cassisi, S., King, I. R., 

Momany, Y., Carraro, G., 2004, ApJ, 605, L125 
Bond, J.R., Arnett W.D., & Carr, B.J., 1984, ApJ, 280, 825 
Carretta, E., et al., 2009a A&A, 505, 117 

Carretta, E.; Bragaglia, A.; Gratton, R.; Lucatello, S., 2009b, 

A&A, 505, 139 
D'Antona, F.; Caloi, V., 2008, MNRAS, 390, 693 
Decressin, T., Meynet, G., Charbonnel, C, Prantzos, N., 

Ekstrom, S., 2007, A& A, 464, 1029 
D'Ercole, A., Vesperini, E., D'Antona, F., McMillan, S. L. W., 

Recchi, S., 2008, MNRAS, 391, 825 
Farrell, S.A., Webb N.A., Barret, D., Godet, O., Rodrigues, J.M., 

2009, Nature 460, 73 
Ferraro, F.R. et al., 2009, Nature, 462, 483 
Gebhardt, K., Rich, R.M., Ho, L.C., 2005, ApJ, 634, 1093 
Glebbeek, E.; Gaburov, E.; de Mink, S. E.; Pols, O. R.; Portegies 

Zwart, S. F., 2009, A&A, 497, 255 
Gratton, R., Sneden, C, Carretta, E., 2004, ARA&A, 42, 385 
Giirkan, M.A., Freitag, M., Rasio, F.A., 2004, ApJ, 604, 632 
Heggie, D.C., Hut, P., Mineshige, S.; Makino, J., Baumgardt, H., 

2007, PASJ, 59, Lll 
Ibata, R. et al., 2009, ApJ, 699, L169 
King, A., 2003, ApJ, 596, L27 
King, A., 2005, ApJ, 635, L121 
King, I., 1962, AJ, 67, 471 

Krolik, J., Active Galactic Nuclei, Princeton University Press 
(1999) 

Kroupa P., Tout, C.A., Gilmore, G., 1993, MNRAS, 262, 545 
Lee J.-W. et al., 2009, Nature, 462, 480 
Lee Y.-W. et al., 1999, Nature, 402, 55 



not straightforward to connect the initial properties of 
a cluster during the growth of the seed black hole with 
the current properties of the cluster hosting an IMBH. 
However, given the loss of most of the FG population 
and the possible near-universality of the seed black hole 
mass, the most likely outcome of the scenario described 
here seems to be one in which the present IMBH mass 
has only a weak dependence on current cluster proper- 
ties. 
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